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Abstract 


Despite the progress that has been made in the measurement of productivity, there have been few 
attempts to cast emissions of bad outputs within a joint-production framework. This paper does 
so. It proposes an experimental framework that depends critically on the shadow price estimates 
of emissions. These shadow price estimates result from the estimation of a cost function for each 
industry of the Canadian business sector. Using a detailed industry dataset that accounts for 
marketed output and greenhouse gases over the 1981-1996 period, our results suggest that the 
shadow values of greenhouse gas emissions are significantly different from zero. Our results also 
indicate that, under a non constant returns to scale technology, failing to account for greenhouse 
gas emissions understates productivity growth. 


Keywords: shadow prices, greenhouse gases, productivity 
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Executive Summary 


Productivity growth is a major source of the improvement in living standards through its positive 
impact on real incomes. However, the link between estimates of productivity growth and 
increases in real living standards cannot be accurately gauged if productivity measures do not 
take into account environmental effects of production. Unfortunately, estimates of productivity 
growth often ignore the environmental impacts of economic activity. 


Ideally, estimates of productivity growth should take account of all inputs and outputs associated 
with a production process, including changes to the environment. In practice, productivity growth 
is normally estimated using techniques that only take account of inputs and outputs that are 
priced. Since most environmental impacts are not traded in markets, they rarely have observable 
prices and so tend to be ignored when estimating productivity growth. 


This paper develops and applies an experimental measure of productivity growth that can 
incorporate environmental impacts. To do so, it estimates a shadow or implicit price of the 
unpriced emissions. The methodology builds on the established technique of a cost-function- 
based model and uses data that relate greenhouse gas emissions and industry output over the 
1981-1996 period. We apply our methodology to one of the more important environmental issues 
facing Canada — greenhouse gas emissions, a relatively novel application. 


We use our methodology to estimate partial abatement elasticities. These show the proportionate 
change in the private cost required to abate greenhouse gas emissions by 1 per cent, everything 
else held constant. Based on past relationships between costs of production and quantities of 
emissions, the estimated elasticities averaged -0.14, with a wide variation across industries. This 
reflects the fact that different industries with different technological structures have reacted 
differently in the past to changes to greenhouse gas emissions. An elasticity of -0.14 indicates 
that a 1 per cent reduction in greenhouse gas emissions requires a 0.14 per cent increase in the 
industry’s private cost. 


Greenhouse gas emissions are found to be substitutes of capital input, indicating that more 
capital is required to reduce greenhouse gas emissions. In contrast, energy input is a complement 
to greenhouse gas emissions, but the relationship between labour input and greenhouse gases is 


Statistically not significant. 


For comparison purposes, we estimated productivity growth both excluding and including 
greenhouse gas emissions in a non-constant returns to scale framework. Overall, our results show 
that ignoring greenhouse gas emissions causes the productivity growth of the Canadian business 
sector to be under-estimated by 21 percent on average over the 1981-1996 period. 


Ee ee cee 
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The framework employed in this paper allows us to take into account the unpriced environmental 
impacts into the productivity growth estimates. The methodology provides useful insights into 
how estimates of productivity growth are affected by the environmental impacts of economic 
activity. A disadvantage is that the methodology is data intensive and technically challenging. 
The productivity growth estimates that the experimental approach produces depends on the 
shadow price of CO emissions. Because it is derived from a multivariate statistical analyses, 
there is an inherent uncertainty in the accuracy of the shadow price. It should be recognized that 
the results of these analyses are therefore subject to error. The size of the error will depend on the 
nature of the functional forms used, the type of econometric analysis employed and the accuracy 
of the data that are utilized. 
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I. Introduction 


Greenhouse gas emissions (GHG: carbon dioxide, methane, nitrous gases and 
chloroflurocarbons) and the prospect of global climate change has spawned a broad research 
program to identify, model, and determine the impact of increases on mean global temperature. 
There is increasing concern regarding the economic and environmental costs of the greenhouse 
gases effect (global warming). Various groups have called for concerted international action to 
reduce emissions of CO. In response to these calls, proposals for emissions cuts have emerged at 
various international conferences, notably the World Climates Conferences at Toronto (1988) 
and Cairo (1990). 


The recognition that economic activity is contributing, through greenhouse gases, to global 
warming, with its attendant costs, has set in motion a substantial research effort. This has taken 
place primarily within the domains of climatological and physical modelling, with widespread 
attention being paid to the economic issues concerning climatic change only in the last decade. 
Early examples of economic analysis were concerned predominantly with the relative efficacy of 
different policy instruments for achieving given targets of emission reduction at the 
macroeconomic level. 


A number of contributions have summarized what we have learned and what we still do not 
know about global warming (Panel on Policy Implications of Greenhouse Warming (1992); Cline 
(1992) and Nordhaus (1994)). While there are many physical, biological, and social dimensions 
to global warming, public policy requires information in at least four areas: a) evidence that 
global warming has occurred or is likely to occur in the future; b) predictions on the magnitude 
and timing of global warming, c) estimates of the damage from global warming, and d) estimates 
of the cost of taking actions to reduce global warming or to mitigate its damage. This paper falls 
in the latter category. Specifically it looks at the productivity effect of reducing GHG emissions 
by treating emissions as part of the production process. 


This paper develops and applies an experimental measure of multifactor productivity growth that 
can incorporate unpriced environmental impacts. This method requires an estimate of the implicit 
or shadow price of emissions: The methodology to estimate this implicit price builds on the 
established technique of a cost-function model of the Canadian business sector and is applied to 
one of the more significant environmental issues facing Canada—greenhouse gas emissions. 


Economists have tried to devise a productivity indicator that takes into account the production of 
bad outputs. Because multifactor productivity measures have been so closely related to changes 
in living standards, it is natural to ask whether net increases (1.e., after taking account of any 
increase in inputs) in marketed outputs are the only things that should ‘count’ as gains to our 
standard living. Pittman (1983) proposed a multilateral productivity index that includes 
undesirable as well as desirable outputs, valuing emissions by shadow prices. This approach 
represents an important step forward in productivity measurement, one that has been followed up 
only recently. His empirical results showed substantial differences between the conventional and 
the proposed extended productivity measures. More recently, the contributions of Fare and 
Grosskopf (1998) and Gollop and Swinand (2001) accept the notion that changes in commodities 
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that are not priced in markets should be considered in evaluating performance and indicate how 
they should be ‘valued’ in the productivity indexes. 


Fare et al. (1993) have attempted to derive plant-specific shadow prices for emissions equal to 
the cost of the desirable output that must be foregone to reduce the undesirable output by one 
unit. A similar methodology has been used by Hailu and Veeman (2000) to construct multifactor 
productivity estimates for the pulp and paper industries. The adjusted multifactor productivity 
growth measures that they report are higher than the conventional. 


In this paper, we develop a basis for integration and extension of the existing literature on the 
impact of bad outputs on economic performance, focusing in particular on their private implicit 
value. The private benefits to producers of using the environment as a free input take the form of 
higher output, or lower input costs for a given amount of production, than if producers reduced 
the “bad” outputs associated with production. That is, lowering emissions implies either 
decreasing marketed outputs (since “bad” outputs are produced jointly with “good” outputs) or 
increasing inputs (from the substitution of energy input or alternative waste disposal). Thus, 
actions requiring reduction of greenhouse gases impose private costs on the business sector. 


Measuring the private shadow values of greenhouse gases and their link to the demand for other 
inputs and other components of the production structure requires a detailed estimable model of 
the production cost. This is accomplished using a rich industry-level (37 industries for 1981- 
1996) panel data set from the Canadian productivity accounts supplemented by data on 
greenhouse gases and natural capital from Statistics Canada’s environment program. 
Experimental estimates on total capital input that combine both produced capital and natural 
capital were also constructed and implemented in our empirical work. Econometric 
implementation of the model with this data base allows estimation and statistical inference of the 
shadow or implicit price of greenhouse gases along with its relationship with input demand, 
technological substitution possibilities and productivity performance. 


We find shadow values of greenhouse gases to be significant, larger for the crude petroleum and 
natural gas industry and increasing in magnitude over time. The results indicate that GHG are 
substitutes of capital input and complements to energy input. But the relationship between GHG 
and labour input is statistically not significant. We also find that, on average, conventional 
estimates of multifactor productivity underestimate productivity performance by about 21 percent 
over the 1981-1996 period. 


The remainder of the paper is organized as follows. Section II is devoted to the theoretical 


framework; Section I gives details on data and econometric implementation; Section IV 
presents the results and their interpretation. Section V provides concluding remarks. 
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II. The Theoretical framework 


1. Costs and Shadow Prices 


Measuring the costs and benefits of greenhouse gases and associated environmental damage 
(emissions for short in the sequel) involves explicitly modeling the production structure, 
recognizing the wide variety of output (revenue) and input (cost) patterns exhibited in the data. 
Our industry-level data set includes information on the production of one “good” output and one 
associated “bad” output (emissions), and the use of four inputs. 


We base our analysis on a cost-function characterization of input demand in the Canadian 
business sector. For empirical implementation, this cost function is augmented by industrial and 
temporal fixed effects to accommodate differences across industries and time periods. This 
detailed modeling framework allows us to explore a rich set of interactions among output 
production, emissions, and input demand. 


More specifically, our total cost function takes the general form G(Y,B,w,D,t) where Y is the 
output; B is the “bad” output; w is a vector of input prices (labour, L; capital, K ; energy, E; 
other intermediate inputs, M ); D is a vector of dummy variables corresponding to fixed effects 
for each industry; and ¢ is a time trend. The total cost function measures the total of the 
expenditures made by an industry for all the inputs it pays explicitly for labour, capital, energy, 
materials and services. If the behavior of the firm indicates that it assigns a negative value to the 
bad output, such implicit cost does not enter the definition of G, although, as discussed further 
below, it may be inferred from it. 


The variable B is included in the cost function on the realization that bad outputs are produced 
jointly with Y , or, conversely, that the environment is used as an input by producers when they 
release effluents into the atmosphere. Production of the bad output allows more marketable or 
good outputs Y to be produced from a given combination of paid inputs, or, alternatively, using 
emissions as an input allows the production of a given amount of Y at lower paid input costs. 


Our approach focuses on private production costs, and emissions are treated either as an input to 
the production process or as an output whose negative value is not fully integrated into firms’ 


revenues. The associated private shadow values zg of the bad output, i.e. the (input) cost saving 


from allowing emissions, may be measured as the cost effects —3¢ = zz. This shadow value 


reflects the marginal amount the producer incurs as a result of a reduction in B. We expect that 
ZB << 0 : 


In our framework, these shadow values incorporate the past behavioral motivations underlying 
cost-efficient production choices, as well as technological substitution possibilities. zz, should 
thus be interpreted as a private cost to producers, since it represents the amount that expenditure 
on other inputs increased (at a given output level) as emissions were reduced by one unit. 
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By Shephard’s lemma, the demand for input j is X; = se (where w, is the market price of input 


j). Then the effect of a change of B on the demand for input j (a second order cost effect) 
measures the dependence of the input use on the ability to dispose of emissions. For example, 


because the bad output is often directly related to the use of energy, decreases in E would be 
associated with declines in B (while increases in most other inputs might be required to reduce | 
emissions). 


The shadow value of emissions: zg = -36 is a relationship analogous to Shephard’s lemma. In 


general, this is a function of all arguments of the G(-) function. It follows from Young’s 
theorem that the impact on zg of a change in an input price is symmetric to the effect of a 
change in B on the demand for X ;: 


OZR 0G 0°G Ox ; 


dw; OBdw, owjdB OB 
2. Scale, Scope, and Various Relevant Technology Properties 


The framework also allows us to gain insights into the existence of cost subadditivity (cost 
savings) that arises from the production of a variety of outputs (economies of scope or diversity) 
or from larger scale of outputs (overall scale economies). Economies of scope, or cost 
complementarities, exist when joint production lowers aggregate costs. Overall scale economies 
(OSE ), a measure of total output cost elasticity as a firm expands outputs along an output ray 


emanating from the origin, holding output mix constant, is defined as 


_ dfnG  dlnG 


ba to EG y eG R- 
Sg RET Ge: 


If OSE is equal to (is greater than) (is less than) 1, production of this vector of output exhibits 
constant (decreasing) (increasing) returns to scale. 


Economies of scope (ESC), a second indicator of the extent to which jointness among outputs 
gives rise to cost savings, can be examined through the second order cost effects 


0’G oe OMCy = dzy 


RCs pay Beige 


The shadow value of the good output Y is its marginal cost 9G = zy ; indeed, if the industry is 
competitive, this should be equal to the price of the good output. Cost complementarities imply 
economies of scope .if ESC <0; diseconomies of scope are indicated by ESC >0, whereas 
ESC = 0 is consistent with the existence of additive costs. 

Neat dLnMC- 
The elasticity €,, 3 =;5° 


provides some indication of private producers’ motivations to adapt output levels to exogenous 
(e.g. regulatory) changes in emissions. For example, under competition, to the extent that MCy is 


measures the impact of emission reduction on marginal cost; it 
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increasing in Y , an increase in MCy resulting from emissions reduction would suggest reduced 
production of Y . Other market configurations can be analyzed as well. 


3. Multifactor Productivity 

The information about the production structure derived in the previous sections is required in 
order to measure and decompose multifactor productivity growth between overall scale 
economies and technical change. As is well known, multifactor productivity growth and technical 


change coincide only under constant returns to scale. 


Total differentiating the total cost function G(Y,B,w,D,t) with respect to time yields 


dG _ yy dw; 0GdY | 0GaB dG a) 
dt 4dw,; dt oY dt OB dt oF 
Dividing “ by G, setting se = X; (Shephard’s lemma) and rearranging the terms, we get 
eS eee 2 
where 
“ = oe = the total change in the cost function, 
(ee ae = the proportional shift in the cost function, 
Egy = Boy. = the cost elasticity of the good output Y , 
0YG 
GB soe = the cost elasticity of the good output B, 


and Q is the aggregate output growth rate , the sum of the good and the bad outputs growth rates 
weighted by their respective cost elasticity, that is 


O = Envy AF EcpB : 


Totally differentiating G= yw ix ; With respect to time and rearranging, yields 


7 arn =< ia X," Substituting this equation into (2) yields 
uF il] 


SG = Ea. ar EnpB a ip (3) 


ee ae ee eee 
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where [| = > ee is the weighted growth rate of all inputs. Following standard practice, we 
os 


define the multifactor productivity growth rate as MEE =Q-I] —that is, the growth of the 
aggregate output not accounted for by the growth of the combined inputs. Substituting, we obtain 
the following decomposition of the multifactor productivity growth rate 

MFP 


eel! Ulsan + Gp )]. (4) 


This formula will be used to compute multifactor productivity and decompose this measure into 
its two basic factors: a) a shift in the cost function due to technical change (-T ); and b) a 


movement along the cost function due to overall scale economies Q[1—(€Gy + Gp NE 


In the above framework, both the good and the bad output are treated as exogenous. Thus Y may 
be set competitively or under imperfect competition; similarly B may be affected by regulations, 
pressures from environment protection groups, etc. 


In order to compute multifactor productivity growth and decompose it into various components 
discussed above, we need data on input and output growth rates, and information on the various 
cost elasticities involved in the above formulas. We can obtain this information from 
econometric estimates of the cost structure G(-). 


III, Econometric Implementation 


1. Specification of the Model 


The cost function implied by the model presented in the previous section takes the general form 
G =G(Y,B,we,w_,We,wy,D;,t), where the general vector representation has been expanded to 
make the individual arguments of the function explicit. The vector of fixed effects D,; includes 
36 industry-specific intercepts with cross effects for each input price and output quantity. 


Econometric implementation of the model and construction of parametric derivative and 
elasticity measures first requires specifying a functional form for G(-). Before doing so we 


express all costs and prices in terms of the price of intermediate inputs, which is a way to impose 
homogeneity in input prices. Then we choose a translog form: 
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G 
a a \: Bon(D;,) + a Bin, (Dy, envin + By, (Dy )enY, + Ban (Dp )£nB, + Bu, (Dit 
+ Byy (EnY, ys + Baan (Dy) (4B, y + De 


+9 > Bin (Dr envinenv jn +) Birn(Dr envinenY, (5) 


izj j i 
oy Pipn( Dp envinlnBy, + Dy Bus (Dy, envint + Byp (Dp )enY,£nB;, 
+ By, LnY,t ate Bp lnBypt. 


The subscripts i and j denote the inputs K,L,and E —capital, labour, and energy—while h is an 
industry index. v;, is the relative input price, defined as v;, = mae where w,,,, is the Fisher price 


index of intermediate inputs.’ Unless explicitly written as functions of industry dummies, as 
discussed below, the #’s are constant parameters. 


Interindustry differences are captured through the following parameterisation of the (’s in (5): 
Bon(Dn) = Bo + ¥,,%orPn » Bin(Dn) = Bi + Y, ,%Dr > Byn(Dn) = By + Y),@mDr » 
Ben(Dy) = Bp + >), 2BnDn » Bin(Dn) = B+ Y), O%nDn + Bin (Dn) = By + Y ,inPr » 

Bin (Dn) = Bir + ¥), Cima » Bizn(Dn) = Bis +), %isnDn and Bin (Dj) = Bir + Y , inDn 


Differentiating (5) with respect to énv,, and using Shephard Lemma gives the share of input 
i= K,Land E in total cost. 


Win = B+ » Bin(Dn )envin + > Bin(Dr )énY;, + > Bixn (Dy, )énB;, + Bir, (Dp dt (6) 


where @;, = Soak . The share of the material and service inputs is calculated as w,,, =1- > ,@in : 


since there are only n—1 independent equations in the model. 


The system of equations (5), along with its associated share equations, (6) should satisfy the 
usual regularity conditions. In particular, for the cost function to be concave in input prices, its 


Hessian matrix | ie | of second-order derivatives with respect to variable input prices should 
Dery Uj 


be negative semidefinite. In addition, the cost function should be nondecreasing in output. Linear 
homogeneity in prices is imposed by construction but can be tested by introducing the price of 
intermediate inputs into each share equation, and adding terms in the cost equation in such a way 
that the augmented share equations can still be obtained from the cost equation by partial 
differentiation as indicated. 


1 The Fisher price indices of material and services were aggregated using the share of materials and services in the 
total cost of these two intermediate inputs. 
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We also assume that the error terms attached to the above equations are optimizing errors and are 
jointly distributed with zero expected value, and with a positive definite symmetric covariance 
matrix. 


2. The Data and Trend Analysis 


A) Methodology 


The model detailed in the previous section is estimated using data for 37 two-digit industries of 
the Canadian business sector during the period 1981 to 1996 (see Table 1 for the list of 
industries).” 


The data set is based on an extended version of the KLEMS database developed on an 
experimental basis to support research themes on productivity and environmental issues. The 
Environmental KLEMS (E-KLEMS) database contains information, at the industry level from 
1981 onwards, on the value of output, the cost of primary and intermediate inputs and Fisher 
chained volume and price indexes of output and inputs. These variables are supplemented with a 
set of environmental variables on waste (greenhouse gases and other pollutants) and natural 
inputs such as natural reserves assets and water use developed by Statistics Canada’s 
Environmental Accounts (Statistics Canada 1997). 


The development of the E-KLEMS database is intended to help address issues on eco-efficiency, 
a measure of the extent to which the economic activity makes use efficiently of the environment 
as a free input. Economic activity has a complex relationship with the environment. It provides 
the raw materials for the production of goods and services that support our living standard, but it 
also causes damage to the environment through the activities of businesses. The conventional 
productivity framework is sometimes criticized for including the value of goods and services 
produced and the income generated through the use of environmental assets, but not reflecting 
the economic cost of depleting those assets or the damage that arises from economic activity. 


There are various features of the E-KLEMS database: 


First, E-KLEMS records the value of environmental assets that are defined as being within the 
scope of the system of national accounts-known as the asset boundary. For an asset to be 
included within the asset boundary of the national accounts it must have an identifiable owner, 
and the owner must be able to derive an economic benefit from the use of the asset. Assets 
included are those termed economic environmental assets such as subsoil assets, land, forests, 
water that are under the control of an economic agent (often the government). 


We only retained the industries for which real output is reasonably accurately measured. Finance and real estate, 
insurance, amusement and recreational service, accommodation and food services, health and social service, 
business service, personal and household service and educational service have been excluded. 
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Table 1. Industries of the Business Sector 


Industries 
] Agricultural and related services 
2 Fishing and tapping 
3 Logging and forestry 
4 Mining 
5 Crude petroleum and natural gas 
6 Quarry and sand pit 
7 Services incidental to mineral extraction 
8 Food 
9 Beverage 
10 Tobacco 
11 Rubber products 
ss Plastic products 
13 Leather and allied products 
14 Primary textile 
15 Textile products 
16 | Clothing 
17 Wood 
18 Furniture and fixture 
19 Paper and allied products 
20 Printing and publishing 
21 Primary metal 
22 Fabricated metal products 
25 Machinery 
24 Transportation equipment 
25 Electrical and electronic equipment 
26 Non-metallic mineral products 
Zh Refined petroleum products 
28 Chemical industries 
29 Other manufacturing 
30 Construction 
31 Transportation industries 
32 Pipeline transport 
33 Storage and warehousing 
34 Communication 
35 Other utility 
36 Wholesale trade 


Retail trade 


Second, the environmental data appended to the conventional KLEMS database are generally 
measured in physical terms. This is the case for greenhouse gases, the bad output retained in the 
present study, and for water use and water discharge exploited in a different study (see Dachraoui 
and Harchaoui, in process). While the value of greenhouse gas emissions is derived implicitly 
using estimates of shadow prices devised econometrically in this study, for the subsoil assets, the 
E-KLEMS makes use of the value of sub soil assets along with the value of the resource rent, 
both of which have been developed by Statistics Canada’s Environmental Accounts. The 
resource rent is the value of the capital services provided by a natural asset. It is calculated as the 
unit value of the output of the natural resource production (e.g. oil and gas) net of the unit cost 
(which includes cost of primary and intermediate inputs). The resource rent in each period is then 
discounted to derive the net present value of the natural asset. 


a rn nnn EET Et en en Le ee nn ee 
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The methodology used in the Canadian Productivity Accounts to construct capital input has been 
amended on an experimental basis to include natural capital as part of the domain of definition of 
capital input. We used the capital stock estimates in current prices for twenty two produced 
assets and two non produced assets, land and natural reserves stocks, along with their respective 


user cost. 


The estimates of the extended capital input are adjusted for compositional changes. In order to 
perform this adjustment, the rental price of these twenty four types of capital are needed. Because 
the rental price is not directly observable, we derive it implicitly on the basis of the available 
information on capital compensation, the value of capital stock assets, capital gains, the tax 
parameters such as the corporate income tax and investment tax credits, the depreciation rate of 
produced capital stock and the depletion rate of natural stock. The internal rate of return, 
calculated residually, ensures the consistency between capital compensation and the cost of 
capital of all produced and nonproduced assets. Over the 1981-1996 period, the effect of changes 
in total capital composition, measured as the difference between the growth rate of total capital 
input and the growth rate of total capital stock, grew at an average rate of 1.1 percent. 


The E-KLEMS database also provides hours worked by industry. Household survey data are used 
to disaggregate total hours into hours worked by different types of workers classified by 
demographic variables such as sex, age, and education. Assuming that workers are paid 
proportionately to the value of their marginal products, Gu et al. (2002) calculate labour input as 
a weighted sum of hours worked by different types of workers, weighted by relative wage rates. 
Annual growth in the labour input for the business sector as a whole from 1981-1996 averaged 
2.1 percent; hours grew an average 1.3 percent per year; and labour composition increased an 
average of 0.8 percent. 


B) Trend Analysis of the Data 


We present certain selective descriptive statistics on the cost and prices of the 37 industries in our 
data. In Table 2, we provide the average levels of total cost and average annual growth rates of 
real output and greenhouse gases, prices of real inputs and cost shares for the period 1981-1996 
for the 37 industries. 


As is clear from the descriptive statistics, the size of the industries, measured by total cost or 
gross nominal output, varies considerably. Trade, construction, transportation equipment, food 
and transportation are among the largest industries in the business sector defined in this paper. 
Other industries such as mining, tobacco, furniture and fixtures, and leather and allied products 
are relatively small. 
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In addition, factor cost shares vary considerably among the 37 industries. For example, labour 
compensations’s share ranges from a low of about 0.09 in crude petroleum and natural gas to a 
high of 0.55 in retail trade. Capital compensation’s share of total cost also varies considerably 
across industries, ranging from 0.09 in construction to 0.67 in crude petroleum and natural gas. 
Generally, capital compensation’s share of total cost, with a few exceptions (most notably 
mining, crude petroleum and natural gas, fishing and trapping, chemical, pipeline transport, other 
utility, beverage and tobacco) is less than labour compensation’s share. Intermediate inputs on 
the other hand, have the largest share in total cost in almost all industries, ranging from 0.12 in 
pipeline transport to 0.90 in refined petroleum and coal products. 


Output grew at 2.5 percent over the 1981-1996 period. The rates of growth of output and inputs 
shown in Table 1 also vary among industries over the period 1981-1996. In leather and allied 
products, tobacco, clothing and refined petroleum products, the growth of output was negative. 
Other industries such as services incidental to mining, fishing and trapping, primary textile, 
textile products, printing and publishing, etc., show a lacklustre growth of output. A number of 
industries experienced output growth rates ranging between 1 percent to over 2.5 percent. Some 
industries in manufacturing and service sectors experienced impressive gains in output; the 
growth rates for these industries ranged from approximately 3.4 percent in rubber products to 
about 6.9 percent in electric and electronic products. The diversity in the growth pattern of output 
and inputs across industries suggests that different industries have experienced different changes 
in their input mix and output and productivity growth. Similar patterns of negative, small and 
rapid growth rates are visible in the growth rates of labour, capital and intermediate inputs. The 
growth rates of output price and input prices with few exceptions were generally positive but 
varied considerably across industries. 


During the 1981-1996 period, the business sector’s greenhouse gases experienced a 1.9 percent 
average annual growth, compared with 2.5 percent for output. Over this period, primary sector 
industries, manufacturing, utility and transportation were the largest producers of greenhouse gas 
emissions. In 1996, the generation of greenhouse gases by these four sectors accounted for 86.8 
percent of the business sector total emissions, unchanged from 1981 (Figure 1). 
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Figure 1. Canadian Business Carbon Dioxide Emissions, 1996 
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Greenhouse gases annual growth rates over the 1981-1996 period also display a great deal of 
variation across industries. About 30 percent of all industries posted a decline in greenhouse gas 
emissions during this period. This includes mostly manufacturing industries such as beverage and 
clothing but also primary industries such as agriculture, mining and services incidental to mining 
which account for reasonably high contribution within the primary sector. Greenhouse gases 
experienced rapid growth in industries such as crude petroleum and natural gas, utility, 
transportation, refined petroleum and coal products. But in the majority of industries, the ratio of 
good output to greenhouse gases grew rapidly, which is a form of efficiency gain. 


The substantial diversity in the growth of output greenhouse gases and the structure of costs 
among the industries over the period 1981-1996 provides a rich body of data to test 
econometrically the impact of different variables on the growth of output and productivity. The 
diversity pattern noted here implies that the response of various industries to changes in variables 
such as greenhouse gases or output are likely to be very diverse. Therefore, we expect the 
estimated elasticities, the shadow price of greenhouse gases and multifactor productivity growth 
rates calculated for different industries using the parameter estimates of our econometric model 
to vary across industries. These inter-industry variations motivate the use of a specification that 


captures industry idiosyncrasies. 
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3. Estimation 


The estimation was carried out for 37 industries of the Canadian business sector for the period 
1981-1996. These industries cover the mineral sector, the manufacturing sector, the non-financial 
services sector and utility industries. In addition to estimates for the 37 industries listed in 
Table 1, we also derived an estimate for these entire sectors by aggregating data across industries 
using Fisher indices. 


The system of equations used to estimate the parameters required by our measurement 
framework consists of the cost function (5); and the share equations for i= K,Land E given by 
(6). The share equation of intermediate inputs is obtained residually from the constraint that cost 
shares must sum to one. We have pooled time-series cross section data for 37 two-digit industries 
of the Canadian business sector for the period 1981-1996 to estimate the model. Estimating the 
model as a pooled system not only adds flexibility to the model (additional degrees of freedom) 
but also imposes cross-equation restrictions to allow a fully integrated cost structure model, 
facilitating more efficient estimates. Seemingly unrelated regressions techniques were used for 
estimation, since the equations share common parameters. 


This equation system was estimated by non-linear seemingly unrelated (SUR) systems 
procedures instead of instrumental variables (IV), which is often used to take into account 
potential output endogeneity or errors in variables. However, we indirectly use the [V technique. 
The instrumental variables technique usually relies on lagged exogeneous variables as 
instruments. Our correction for first-order autoregressive disturbances makes use of lagged 
values of exogenous variables as instruments. Therefore, the correction for autocorrelation of 
residuals turns out to rely implicitly on the IV technique. 


Adaptations were made to accommodate potential unknown sources of heteroskedasticity. 
Changing the input demand equations to input/output measures to reduce variations in scale 
across industries and time did not affect the estimates substantially. The results reported below 
are based on the unmodified system, using White’s heteroskedastic-consistent covariance matrix 
to generate standard errors. 


Durbin-Watson tests indicated that autocorrelated errors were present in the cost and input 
demand equations. Therefore, the lagged dependent variable was incorporated into the cost 
equation giving the form G,=a@+fX,+G,,+u, which in tur implies, after suitable 
substitutions, the following form G, = a(1+p)+ B(X, + PX, )+(u, + pu, ) where X, refers to 
the vector of right-hand variables in (5), @ and f refer to the corresponding parameters, and p— is 
the coefficient of autocorrelation. 
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IV. Empirical Results 


The estimated parameters and their standard errors as well as the parameter estimates for industry 
dummies are shown in Table 3. Although in a model this complex, the individual parameter 
estimates have limited interpretation, the overall statistical significance of the parameters is 
notable. Most industry dummies are significant; so is the estimate of 9. The R*s, reported in 
Table 3, indicate excellent “fits” for the estimated equations—all being higher than 0.95. The 
results also indicate that the model is well estimated. The square of the correlation coefficients 
between the actual and predicted values is high, and the standard errors of each equation are 
small. In addition, all the required regularity conditions are satisfied at each point in the sample. 
The coefficients of the model are statistically significant and have the correct sign. 


I. Specification and Hypothesis Tests 


The results of the hypothesis tests using log-likelihood ratios are shown in Table 4. The results of 
the hypothesis tests using log-likelihood ratios decisively reject the joint hypothesis that the 
dummy industry coefficients are zero (first row), indicating that strong interindustry differences 
are present in the cost structure of the industries under consideration. The hypothesis that the 
coefficients of the bad output are zero in the total cost function (5) is decisively rejected (second 
row). Also, the hypothesis that the cost function is homogeneous of degree one in the input prices 
is not rejected. 


Similarly, the hypothesis that firms do not produce a joint product é¢, =0, that there is no 
technical change prevail, €¢, =0, and the industries operate under overall constant returns to 
scale Egy + €gg =1 were separately tested. In each case, the test consists in obtaining a vector of 


estimated parameters © from the four equation system (5)-(6) plus the additional restriction to be 
tested. If © is the vector of unrestricted estimated parameters, then the quadratic form 


Me= (9 — 6) {Cov(®) _ Cov(6)} "(6 = 9), 


is asymptotically chi squared with degrees of freedom equal to the number parameters of the 
imposed restriction. Each of the first two hypotheses was rejected, but the assumption of constant 


ray scale economies was not rejected (M =210> 765.09) =134;M =157 > 745.001 = 134, 


M =108 > 765.001 = 114 respectively). 


oe en De EE eee 
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Table 3. Parameter Estimates of Cost Structure (Sample Period 1981-1996) 


Parameter Estimate Standard Error Parameter Estimate Standard Error 
Bo 1.79844 2.22285 OK 14 -0.30771 0.241422 
ao aK15 -0.3605 0.204562 
ao 1.68764 1.71969 aK. 16 -0.35499 0.200541 
G0 3 1.19111 1.63142 AK 17 -0.39906 0.1789 
G04 -0.03073 1.65079 OK 18 -0.14448 0.197741 
to 5 4.09267 1.66929 ax 19 -0.13982 0.182936 
ao7 0.942586 1.69085 OK2I -0.38251 0.200502 
aos 1.76937 2.4728 digas -0.23782 0.196125 
Go.9 22212 2.12094 OK23 -0.28218 0.185983 
Qo.10 1.29155 1.73653 OK 24 -0.48356 0.183295 
Qo11 2.6929 1.69556 OK25 -0.18819 0.218136 
Go.12 2.06508 1.58799 aK.26 -0.17009 0.178373 
a3 -1.38548 1.7953 Gest -0.20937 0.188203 
Go 14 1.75288 1.92197 aK 28 -0.48084 0.236293 
9,16 2.22914 1.77913  AK30 -0.23059 0.183778 
Qo.17 1.1383 1.62589 AK, 31 -0.60302 0.282178 
Qo 18 0.976801 1.62156 ards 0.129868 0.201813 
0,19 0.183455 1.63877 OK 33 1.24E-03 0.223328 
0,20 1.08851 1.7916 OK 34 -0.55209 0.229791 
0,21 1.93164 1.77657 AK 35 0.068963 0.478113 
Q0,22 1270147, 1.67682 OK 36 -0.44587 0.213077 
Q0,23 0.8495 1.64661 OK 37 -0.40323 0.193961 
10,24 3.27164 1.63606 
Q0,25 2.26088 1.6426 Br 0.36534 0.238352 
Q0,26 -0.03795 1.69182 GL] 

00,27 1.76675 2.14642 aL2 0.37586 0.238447 
0,28 2.74055 1.92493 QL,3 0.847295 0.281046 
010,29 1.02683 1.60697 a4 -0.06589 0.324054 
20,30 0.68629 1.67595 GL,s -0.19844 0.287238 
Qo,31 3.23662 1.71144 aL6 -0.06644 0.266611 
0,32 4.97099 1.73975 a17 -0.06163 0.279595 
0,33 1.41906 1.87904 ang 0.219847 0.392814 
O09, 34 4.92566 1.62951 az,.9 0.377695 0.275679 
Qo,35 -0.1748 wievioy GL, 10 0.419143 0.256785 
0,36 2.57108 1.65373 QL11 0.165805 0.269847 
Qo, 37 2.68888 1.66999 Gr, 12 0.254067 0.280473 

OL, 13 0.323133 0.291701 
Br -0.13792 0.188758 QL, 14 -0.09217 0.243175 
a | QL15 0.033168 0.254252 
ax2 -0.48223 0.175482 GL, 16 0.182877 0.332996 
OK3 -0.21681 0.172512 OL,17 0.116462 0.309667 
aK4 -0.29212 0.178739 GL, 18 0.152661 0.27601 
aks 0.524024 0.189992 ; OL,19 -0.40298 0.257467 
axe -0.1806 0.168562 1,20 0.10033 0.290089 
ax7 -0.29378 0.194858 G12] -0.05878 0.246533 
ans -0.3274 0.463629 Q1,22 0.012101 0.266959 
ano 0.013891 0.225846 G1,23 0.238059 0.251279 
AK, 10 -0.57332 0.201513 G124 0.211646 0.285517 
OK 11 -0.34389 0.184712 OL25 0.232119 0.271288 
AK, 12 -0.23037 0.196765 1,26 -0.01764 0.270478 
OK 13 -0.03914 0.2488 O1,27 -0.04786 0.239775 
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Table 3. Parameter Estimates of Cost Structure (Sample Period 1981-1996) - Continued 


»*Parameter 


OL, 28 
Q1,29 
aL, 30 
QL 31 
GL, 32 
GL, 33 
QL 34 
OL, 35 
aL 36 
OL, 37 


BE 
QE 1 
OE2 
aE 3 
OE, 4 
GE s5 
aE 6 
On7 
aE8 
QaE9 
QE, 10 
QE 11 
OE, 12 
GE 13 
GE, 14 
OE 15 
GE, 16 
OE 17 
QE, 18 
QE 19 
QGE,20 
QF 21 
OF, 22 
GE 23 
QE 24 
GE 25 
GRE 26 
QE27 
GE 28 
QE, 29 
OE 30 
OF 3] 
OE, 32 
GE, 33 
OR 34 
OE 35 
GE 36 
GE 37 
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Estimate 
-0.14453 
0.206469 
0.170548 
0.067362 
0.378274 
0.154204 
0.549666 
0.516156 
-0.50552 

0.54435 


0.076285 


0.021359 
-0.01286 
0.01528 
-0.03556 
0.026359 
-0.01064 
-0.04435 
-0.0402 
-0.04314 
-0.03571 
-0.03449 
-0.05312 
-0.03205 
-0.04082 
-0.05191 
-0.02588 
-0.04681 
0.013567 
-0.04913 
0.019871 
-0.04401 
-0.05048 
-0.04214 
-0.05461 
1.14E-03 
-0.02463 
9.53E-03 
-0.04953 
-0.04157 
0.0548 
-6.36E-03 
-0.0176 
-0.03765 
0.017356 
-0.01819 
-0.01638 


-1.9234 


-0.1233 
-2.61E-03 


Standard Error 


0.318454 
0.266481 
0.277321 
0.365111 
0.265955 
0.282896 
0.361923 
0.406778 
0.439498 
0.363137 


0.018176 


7.72E-03 
7.45E-03 
7.57E-03 
7.67E-03 
7.53E-03 
7.80E-03 
7.63E-03 
7.52E-03 
7.86E-03 


:7.69E-03 


8.01E-03 
7.73E-03 
7.56E-03 
7.56E-03 
7.66E-03 
7.86E-03 
7.72E-03 
7.92E-03 
7.90E-03 
7.65E-03 
7.52E-03 
7.63E-03 
7.84E-03 
8.22E-03 
7.47E-03 
7.68E-03 
7.74E-03 
7.84E-03 
7.86E-03 
7.81E-03 
8.42E-03 
7.68E-03 
8.00E-03 
8.12E-03 
7.93E-03 
7.87E-03 


0.698464 


0.202763 
0.203294 


yh 


Parameter 
Oy 4 
ay 5 
aye 
ay7 
ays 
ayo 
Gy 10 
Oy 11 
Gy 12 
Gy 13 
Qy, 14 
Qay 15 
Gy 16 
ay 17 
ay 18 
ay 19 
Qy,20 
ay, 21 
ay, 22 
ay 23 
Ay 24 
Gy, 25 
Gy 26 
Gy, 27 
ay, 28 
ay, 29 
ay 30 
Oy, 31 
Oy 32 
ay 33 
ay 34 
Ay 35 
Ay, 36 
Ay, 37 


Be 
- OBl 
aR 2 
aB3 
OB4 
OBs5 
OB 6 
OB,7 
oR8 
aBR9 
QB, 10 
OB 11 
QB 12 
OB, 13 


Estimate 
-0.42131 
0.577569 
-0.08828 
0.071283 
-0.01039 
-0.2338 
0.325122 
-0.23438 
-0.19271 
0.336623 
-0.04012 
0.041399 
-0.14175 
-0.02258 
-0.05748 
0.212148 
-0.2662 
-0.03067 
-0.1445 
0.034422 
-0.21407 
-0.33133 
0.181484 
-0.11886 
0.033681 
0.06145 
0.152294 
-0.02587 
-1.30088 
-0.28925 
-0.5157 
-0.18175 
9.02E-03 
-0.24469 


-0.24357 


0.22919 
-0.04362 
0.71505 
-1.92779 
-0.12439 
0.031471 
-0.0389 
-0.03232 
-0.02241 
-1.41E-03 
-0.01937 
5.48E-03 
-0.02948 
-1.73E-03 
0.021641 
0.168601 
1.46E-03 


Standard Error 


0.225944 
0.215312 
0.188994 
0.191973 
0.313918 
0.270028 
0.2155 
0.189483 
0.197137 
0.245626 
0.206837 
0.209288 
0.23685 
0.21998 
0.201745 
0.203255 
0.302603 
0.202297 
0.214563 
0.187276 
0.208113 
0.21013 
0.209717 
0.375543 
0.236817 
0.210291 
0.211514 
0.327246 
0.250253 
0.242029 
0.203595 
0.28968 
0.267514 
0.232764 


0.349497 


0.091914 
0.082894 
0.156504 
0.126386 
0.081974 
0.089945 
0.150138 

0.13163 
0.099737 
0.089924 
0.109947 
0.106752 
0.087042 
0.086581 
0.098459 
0.101881 
0.092735 
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Table 3. Parameter Estimates of Cost Structure (Sample Period 1981-1996) - Continued 


Parameter Estimate Standard Error Parameter Estimate Standard Error 
apie -0.10597 0.111986 B, 0.042821 0.017031 
Gp 20 0.040838 0.140766 Bx -0.02496 7.10E-03 
ap2] 9.37E-03 0.172539 Bur 0.033634 3.96E-03 
Op. 22 0.018946 0.11355 Bee 0.018073 1.26E-03 
Op 23 0.070726 0.107867 Byy 0.233183 0.087873 
OB, 24 2.92E-03 0.131835 Bop 6.74E-03 0.032173 
Ap 25 0.023261 0.090364 Ba -2.94E-04 1.46E-04 
AB 26 0.015604 0.15398 
Gp27 -0.05369 0.141198 Bex 6.24E-03 2.71E-03 
anos -0.1261 0.152542 Bex -7.68E-03 3.13E-03 
gee -9.45E-03 0.085635 bey -0.0164 2.88E-03 
Ap, 30 -0.07449 0.128129 Bes 4.47E-03 1.39E-03 
OB 31 -0.0695 0.235423 Ber 1.58E-04 2.23E-04 
Op, 32 0.106971 0.103081 Bri 0.031343 6.91E-03 
i's -0.02328 0.084992 Bry 0.162921 0.014872 
OB, 34 0.016713 0.095885 Bre -0.01553 7.15E-03 
OB 35 0.156925 0.138511 Br -2.83E-03 1.24E-03 
OB 36 -0.11913 0.144599 Bry -0.11251 7.49E-03 
OB 37 0.071887 0.140659 Bip 1.39E-03 3.61E-03 
Bu -5.61E-04 5.91E-04 
fa) 0.2096 0.040871 Byp 0.036789 0.083814 
De 0.810904 0.024865 By -5.14E-03 4.77E-03 
Ox 1.00025 4.02E-03 Ba 3.09E-03 2.99E-03 
01 1.00166 1.46E-03 
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(Sample Period 1996-1998)—Continued 


Table 3. Parameter Estimates of Industries’ Cost Structure 


Equation Standard Error R2 
Total Cost 0.027976 0.98 
Capital Share 0.025634 0.98 
Labour Share 0.012446 0.99 
Energy Share 0.004354 0.98 


Table 4. Hypothesis Testing 


Parameter Restrictions Log of xe Degrees of 
Likelihood Freedom 
Cho Skies eLh a Rhee. Sb See ee ety hoy 194] 1239 407 
Qn, = 2 = 9 
1863 cay. 64 
2B h Qi Bh = 0 


Note: The critical values x? with 407 and 64 degrees of freedom are 654 and 126, respectively. 


q@ is the vector of dummy parameters. 
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2. Various Estimation Issues 


A) Spurious Correlation 


The presence of common trend among variables in the time series models of production structure 
is a serious econometric issue. This criticism is equally applicable to production and cost 
function studies, whether they include environmental variables or not. It is true that private sector 
variables such as output, labour, intermediate inputs and private capital stock are highly 
correlated over time and may share a common trend. There is nothing particularly different about 
environmental variables in this respect. 


One method for removing a common trend is to estimate the model in a first-difference form. 
Estimation of this form eliminates a potential influence of trend. Equations (5) and (6) were 
estimated in ‘first-difference’ form by setting the serial correlation parameter p to unity. The 
parameter estimates (not reported here, but available on request) indicate that the models fit the 
data very well. Signs and magnitudes are similar to those when the models were estimated in 
level form. This should not come as a surprise as the values of the serial correlation coefficients 
p Shown in Table 3 are close to unity. 


B) The Issue of Capital Fixity 


The various cost and demand relationships developed above are characterized through first and 
second order derivatives or elasticities of the cost function with respect to the arguments of G(-). 
However, divergence in input demand patterns from those appropriately represented by 
Shephard’s lemma would complicate or preclude the estimation and interpretation of these 
elasticities. Even though such deviations from standard assumptions of basic microeconomic 
theory did not emerge in the end, the knowledge that they might exist stimulated an empirical 
investigation of alternative models that recognizes these potential difficulties. 


The most common problem of this sort is the quasi-fixity of inputs such as capital. If full 
adjustment to equilibrium input levels does not take place within the time frame of the data, 
Shephard’s lemma will not appropriately represent input demand behavior. This rigidity problem 
is often dealt with by incorporating capital stock instead of its rental price in the G(-) function if 
we have reasons to believe that capital has binding fixity constraints. This implied divergence 
from equilibrium demand (or, equivalently, variations from full utilization) represented by the 
deviations between a factor’s shadow value zx = -2 and its market price wz. 


Alternatively, the true/effective quantity demand of an input may be represented by directly 
adapting the data to embody the discrepancy. In particular, if the true (or shadow) price of the 
factor zy is used as an argument of G(-)rather than an unadjusted market price, the validity of 


Shephard’s lemma is maintained. 
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Although the data for this study were carefully constructed to reflect the input flow values, 
sensitivity checks were carried out to determine the validity of the assumption of variable inputs. 
These checks supported our final empirical specification; the assumption seemed justified by the 
appropriate levels and shapes of the resulting demand equations. In fact, when capital was not 
characterized as a choice variable, the results were not as justifiable as when Shephard’s lemma 
was implemented. 


Our empirical findings based on these data suggest that the approach employed was carried out in 
a manner consistent with economic theory. The use of Shephard’s lemma seems justified by both 
the correct (in terms of required regularity conditions) and intuitively plausible estimates of 
demand behavior. And when optimization equations were not imposed for the K input, the 
resulting estimates remained substantively unchanged. 


3. Economic Interpretation 


The various cost elasticities computed from the estimated parameters for the full data sample are 
presented in Table 5. The reported estimates are weighted averages across all the 37 industries 
and time periods for each measure. The t-statistics are based on computation of the measures 
evaluated at the average (mean) values of the data.” 


The primary measures for evaluating the marginal benefit of using the environment for disposal 
of greenhouse gases is the shadow value sg computed by partial differentiation of (5) with 
respect to £nB. At the aggregate level, the cost elasticity of the bad output €,, is -0.14, thus 
indicating that allowing higher emissions is cost-saving for the producer. The significance level 
is better than 2 percent. There is however a great deal of variation across industries. With a value 
of -1.91, the oil and gas industry shows the highest cost elasticity with respect to emissions, 
followed far behind by chemical industries with -0.12. 


The negative aggregate estimate of €x 3 suggests that capital has a tendency to “substitute” for 


environmental quality, in the sense that additional capital is required to reduce emissions. Both of 
the elasticities of labour and energy with respect to emissions are positive, but only the latter is 
statistically significant. Reduction of emissions implies lower energy use. 


For the outputs, the positive elasticity (€yc,,; = 0.037 ) implies the absence of scope economies or 


jointness between the good and bad output. However, this result is not statistically significant. 
Some interesting implications are suggested by the emission shadow value in Table 5. The results 
indicate that the shadow value of emissions increases at 0.03 percent per year on average, holding 
other variables constant. 


> The measures were constructed for these data using the delta method (essentially a generalized Wald test) by the 


ANALYZ command in TSP. 
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It is also informative to compare some of the elasticity estimates across time periods and 
industries. We report measures separated into pre- and post-1990s, in an attempt to identify a 
possible break in public awareness or institutional pressure with regard to the effects of 
greenhouse gases. If anything, the €,, value for the pre- and post-1990s show a slight reduction 


in the proportional cost savings of B disposal: €, , = —0.1504;— 0.1208 , respectively. 


The estimates of multifactor productivity growth rates require estimates of ray scale economies, 
aggregated good and bad outputs and technical change. The required cost elasticities are taken 
from the parameter estimates of the translog cost function discussed earlier. Figure 2 presents the 
multifactor productivity indices for the period 1981-1996 for the 37 industries considered in our 
empirical analysis. The estimates used in figure 1 are weighted averages of all industries in the 
sample. In other words, multifactor productivity indices for each industry were estimated and the 
weighted average of these estimates is presented in Figure 2. The standard framework that 
excludes bad output underestimates productivity growth by almost 0.5 percentage points per year 
over the 1981-1996 period.* 


Figure 2: Multifactor Productivity Compounded Growth Rate Under Alternate Production 
Frameworks: Business Sector 


150.00 


Alternative: Goods and bads are jointly included in the technology; 
Conventional: Only goods are considered in the technology. 


140.00 


130.00 


90.00 
1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 


eum A|ternative Framework === Conventional Framework 


Under the constant returns to scale assumption that is consistent with the growth accounting framework, the 
exclusion of greenhouse gases emissions leads to an under estimation of 0.2 percentage points (see Harchaoui et 
al 2002). 
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V. Concluding Remarks 


To a remarkable extent, environmental protection is generally perceived in the public debate as 
imposing costly burdens on the economy, stifling innovation and lowering productivity. 
However, the conclusion that environmental protection generally leads to lower productivity 
performance is in fact an artifact in the way the productivity measure is implemented—a 
methodology that counts only the cost of environmental protection but ignores the production of 
a better environment, say in terms of emissions reduction. 


For many years, many studies hampered by this methodological shortcoming, concluded that in 
general environmental protection leads to a decline in the productivity performance. However, a 
recent strand of the economic literature, to which this paper belongs, recognizes that some 
outputs are valuable when sold and others are damaging when released (see Murtough et al. 2001 
for a review of this literature). When the productivity framework considers an industrial process 
in its entirety, environmental protection is no longer seen to necessarily hamper productivity 
performance. 


This study uses a detailed model of the production structure in the Canadian business sector to 
measure the private costs that producers have incurred in the past two decades as they reduced 
greenhouse gas emissions. We find the private implicit or shadow value of emissions to be 
significant, larger for the mining sector and increasing in magnitude over time. Firms do not 
choose their production and output mix as if the environment was free and valueless. This means 
that firms incur, or perceive, costs from GHG emissions, beyond their private input bill. These 
costs may have resulted from regulatory pressure or from public opinion, or from anticipations of 
future regulatory or public pressure. 


Failing to account for greenhouse gas emissions generally leads to an underestimate of 
productivity growth. Cost changes associated with emission reductions are interpreted as 
productivity losses in conventional measures. We have computed a multifactor productivity 
index that includes emissions as input. That index grows faster by half a percentage point a year 
over the 1981-1996 period than the conventional index. 


In conclusion, we have shown that the methodology used in this paper has the major advantage 
that it can readily incorporate unpriced environmental impacts into productivity growth 
estimates. A disadvantage is that the methodology is data intensive and technically challenging. 
The value of the new productivity estimate that the experimental approach produces depends on 
the shadow price of CO2 emissions. Because it is derived from a multivariate statistical analyses, 
there is an inherent uncertainty in the accuracy of the shadow price. It should be recognized that 
the results of these analyses are therefore subject to error. The size of the error will depend on the 
accuracy of the functional forms used, the type of econometric analysis employed and the 
accuracy of the data that are utilized. 


Nevertheless, the approach utilized in this study can provide useful insights into how estimated 
productivity growth can be affected by the environmental impacts of economic activity. There 
may be scope to extend our analysis to incorporate other environmental by-products. This would 
require, for example, the measurement and consideration of other pollutants. 
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